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Abstract 

We introduce a new fundamental property of waveguides induced by the forces of the guided 
light, namely, the ability to self align or be in instability. A nanoscale waveguide broken by an 
offset and a gap may tend to self align to form a continuous waveguide. Conversely, depending on 
the geometry and light polarization, the two parts of the waveguide may be deflected away from 
each other, thus being in an unstable state. These effects are unique as they rely on the presence 
of both the guided mode and the scattered light. Strong self alignment forces may be facilitated 
by near field interaction with polarization surface charges. 
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I. INTRODUCTION 



Laser light has significant mechanical effects on microscopic objects, as was initially 
pointed out about four decades^ago In addition to the vast work on trapping and 



manipulation of small particles [2|, much effort has been directed at cavity based optome- 
chanical devices, where optical forces may be considerably enhanced 4, sl- A less explored 
option, however, is optical forces on waveguides, which become observable at the microscopic 
scale p, lZII . 

Research of optical forces on waveguides is motivated by the growing capabilities of 
nanofabrication that enable new possibilities of nanoscale light manipulation 

lio|- 



ude the laws of the mechanical 
, [l5, Q. These new physical 
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The theory of guided light is thus being extended to inc 
effects of light on the guiding structure itself [ll, 12, 
mechanisms relate the properties of the guided modes to the forces created by them. For 
instance, light guided between two waveguides or mirrors creates a repulsive force for an 



antisymmetric transverse field and an attractive force for a symmetric transverse field 12, 



Moreover a .uperposition of a,_.c and a„ti.y..etnc ,™d. .ay hold the 



waveguide in a stable equilibrium [13|, ll5|, [ITJ. Such phenomena can be experimentally 



observed in nanomechanical devices fabricated on a chip, as was recently demonstrated with 

n 

a suspended silicon waveguide |18l ]. 

In this study, we introduce a new fundamental property of waveguides, namely the ability 
of a waveguide to self align by light forces when it is perturbed by a small offset misalignment. 
We show that the size of the waveguide and the type of eigenmode determine whether the two 
misaligned parts will tend to self align or deflect away from each other. This phenomenon 
is unique as it relies on both the guided waveguide eigenmode and the scattered radiation 
from the perturbation. A strong self alignment force is created due to polarization surface 
charges that dominate the optical force with a near field interaction. We further investigate 
the exerted forces when a gap is introduced between the two parts of the waveguide. 

The geometry under consideration is shown in Fig. [H A single mode slab waveguide 
of half-width d and permittivity e,- is broken by an offset in the x axis. A, and a gap in 
the z axis g. An eigenmode is incident from the left (input waveguide) carrying power Pin, 
most of which is transmitted to the output waveguide (Pout)? while the remainder is either 
scattered (Psca) or reflected back into the input waveguide. No variations of the geometry 
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FIG. 1: A slab waveguide broken by an offset and a gap. 
are assumed along the y axis, and therefore all quantities given are per unit length. 

II. A WAVEGUIDE PERTURBED BY AN OFFSET 

First, we examine the case of no gap, (7 = 0, and a small offset A. For the calculation of 
the fields we take an approach of mode matching approximation, similar to that described 
by Marcuse Q, [2^. We begin by considering a transverse electric (TE) incident mode, for 
which the nonzero field components are Ey, H^, and Hz- The even guided mode incident 
from the left is given inside the dielectric slab by E[ = Ate cos(A;^x) exp(— j/?gz), where 



Ate = \/'2^fJ'oPm/l3g{d + 7^-*^), is the transverse wavenumber in the dielectric, /5g is 
the longitudinal wavenumber, 7 is the transverse decay constant outside the slab, and the 
time dependence is of the form e-^^*. The transverse electric field in each region may be 
represented as a sum of the guided modes and the continuous spectrum of the radiation 
modes, corresponding to the scattered light. Explicitly, assuming the interface plane is at 
z = 0, then for z < (input waveguide) this field is given by 

/•oo poo 

Ei = E; + a,E, + / dpgci(p)£'ci + / dpgoi(p)£'oi , (1) 
JO Jo 

whereas for ^ > (output waveguide) it reads 

poo poo 

E2 = atEt+ / dpge2(p)^c2+ / dpgo2(p)^o2. (2) 
Jo Jo 

In the above two equations, Or and at are the reflection and transmission coefficients of 
the guided mode, respectively; E,- and Et are the reflected and transmitted guided modes, 
respectively; gei,2 and goi,2 are the amplitudes of the even and odd radiation modes, respec- 
tively; p is the transverse wavenumber of the radiation modes outside the slabs; -E'ei,2 and 
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i?oi,2 are the even and odd radiation modes, respectively. The amphtudes of the guided 
and the radiation modes may be approximated analytically by expressions containing the 
overlap integral between the respective mode and the incident guided mode. 

Generally, force densities on dielectrics may be viewed as resulting from two pro- 



cesses 
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22|: (1) the interaction of effective polarization volume current densities 
with the magnetic field, and (2) the interaction of polarization surface charge densities 
with the electric field. In the TE case, no component of the electric field is perpendicular 
to the boundaries of the dielectrics, and therefore no polarization surface charge densities 
are formed. Thus the time-averaged volume force density in the x direction is given by 
^Re (jcuPy X fiQHz*), where Py = £:o(^r ~ ^)Ey is the polarization density. By integration 
over the volume, we obtain the total transverse force on the output waveguide in terms of 
integration over only the top and bottom surfaces 

= l^oie. - 1) Hdz {\E,\\^^^^ - \E,\\_^^^) . (3) 

When there is no offset (A = 0), the field is symmetric around x = and the guided mode is 
not disturbed by a discontinuity, and thus no force acts on each of the waveguide parts. Once 
an offset is introduced, symmetry is broken, and scattering occurs from the discontinuity. 
The question then arises whether the two parts of the waveguide will tend to deflect away 
from each other, or self-align to form back a continuous waveguide while maximizing the 
output power. 

For the evaluation of the force on the output waveguide, the field expression of Eq. ([2]) 
is substituted into Eq. ([3]). At this point, we are interested in small offsets, so that only 
terms up to the first order of A are kept. Noting that cross-products of the even radiation 
modes with the odd radiation modes are of order larger than A, and using the symmetry 
properties of the modes, the expression for the force reads 



x2 



poo poo 

eo{er-l)Re dz E^* dpgo2(p)^o2, (4) 
Jo Jo 



where the integration is performed a.t x = d + A. Hence, it is evident that the odd radiation 
modes created by the scattering are responsible for the transverse force on the waveguide. A 
direct measure of the waveguide's tendency to move either way is the derivative of the force 
with respect to A, dFa.2/dA(A = 0) denoted by -Fx2'- Bearing in mind that -Fx2(A = 0) = 0, 
the force may be approximated by Fx2^- The only term in the above equation that depends 



on A is qo2i and an analytic expression for dgo2/dA may be obtained. The integration over 
z is then performed analytically, and the closed form expression for the derivative reads 

Jo 1 + ^ cot^ (ad) T + /?r(/5r - /5g) 

where ko = u/c. While the spectrum of radiation modes contains both propagating modes 
having real and evanescent modes having imaginary /3j-, the above expression shows that 
it is only the evanescent radiation modes that contribute to the generation of this force, i.e., 
the integrand is nonzero only for p > /cq. 

Both the force [Eq. (jl])] and the quantity -Fx2'A [Eq. ([5])] are plotted in Fig. [2]^a) as a 
function of d, for an offset of 2% of d. The range of d shown is 20 nm to 110 nm, where the 
slab is single mode in each polarization. The wavelength is taken to be A = 1.55 fim and the 
permittivity is 3.48^, corresponding to silicon at that wavelength. The forces are normalized 
by Fq = Pin/ c, which is the momentum per unit time carried by a plane wave. These results 
are compared with a Finite Element Method (FEM) simulation, and as seen the three 
curves are virtually identical. In addition to the integration of the force on the polarization 
densities, we have also integrated over the Maxwell stress tensor 22I to obtain the force, and 
found excellent agreement between the two methods, which are mathematically equivalent 
for exact solutions of Maxwell's equations. 

The plot of Fig. [2](a) reveals two regimes; the first is a self- alignment regime up to a slab 
half-width of about 57 nm, for which Fx2' < corresponding to a restoring force. The second 
is an instability regime in which a small offset results in a deflection force. Both regimes 
exhibit an optimal slab width for which the force is strongest. Although we show here only 
the transverse force on the output waveguide, when the offset is small, the scattered power 
is negligible and by virtue of momentum conservation, the force on the input waveguide is of 
the same magnitude and opposite in sign to that on the output waveguide. In fact, we have 
shown analytically that F!^^=-F^2- To better illustrate the different regimes, we depict in 
Fig. [3](a) the transverse force Fx2 as function of A for d = 39 nm where a negative restoring 
force is seen, d = 57 nm where the force derivative at A = vanishes at the transition 
between the two regimes, and for d = 110 nm where instability in the form of a deflecting 
force is observed. 

When the incident mode is transverse magnetic (TM), the situation is considerably more 
involved. The field components for the TM mode are E^, E^, and Hy, and the incident 
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FIG. 2: Transverse force as a function of d for an offset A of 2% of d. The permittivity of the 
waveguides is = 3.48^, and the wavelength is A = 1.55 /im. (a) TE incident mode, (b) TM 
incident mode. 




FIG. 3: Transverse force ¥^2 as a function of the offset A for different values of waveguide half- 
width d. The solid line corresponds to the analytic analysis and the square markers indicate FEM 
simulation results, (a) TE incident mode, (b) TM incident mode. 

magnetic field is given by = A^m cos^k^x) exp(— j/5gz), where 

The fields in each region are described by Eqs. ([1]) and ([2]) with E replaced by H. The force 
mechanism in the TM case differs substantially than that of the TE case, as for the TM 
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the electric field has normal components to discontinuities in the dielectric, and therefore 
polarization surface charge densities are formed. These surface charge densities interact with 
to give rise to surface force densities. Considering again the output waveguide, the first 
contribution is from the polarization surface charge densities created by -E^ at the z = 
interface. The force density of this term is given by |Re[— £:o(£^r ~ ^)EzEx*], where is the 
field just inside the slab at 2 = 0. Calculation of this force on the output waveguide while 
keeping only terms up to the first order of A yields 

/ {KPt + c^Pg) sin[(fc^ - (r)d] ^ {k^(3, - (j(3^) sin[(fc^ + a)d] 

where an analytic expression for go2 may be obtained. A second contribution is from the 
polarization surface charge densities formed on the top (x = d+A) and bottom (x = —d+A) 
parts of the waveguide, given by the plus and minus of ^£o(^? ~ respectively. In both 

cases Ex is the field just inside the slab. Summing up the two contributions and integrating 
over z results in the expression 

Fx2,UD ^ ATMpg-^ — cos{k^d) Re / dpgo2(p) x sin(crrf)/ [j{p, - /3g)] . (7) 
[Loeoej.) Jo 

The third contribution comes from the volume force density given by |Re[— eo(^r ~ 
l)jujEzfJ,oHy*], and integration over x and z gives 

Fx2,v ^ -Atm^^^^ — — cos{k^d) Re [ dpgo2(p) x sm{(Td)/ [j{f3, - /3g)] . (8) 

The total force is the sum of all three contributions, and its derivative is obtained by 
analytically differentiating qo2- Similarly to the TE mode, only the evanescent part of the 
odd radiation mode spectrum participates in the generation of the force, and the relation 
F^^ = holds as well. 

The quantities Fx2 A^ Fx2, as well as the force calculated by FEM, are shown in Fig. [2]^b) 
as a function of d for an offset A of 2%. Contrary to the TE case, here instability occurs 
for small values of d, and above about 70 nm there is self-alignment. The restoring force 
increases monotonically, so that at the maximum value in the shown range of d, it is about 
two orders of magnitude stronger than the peak value of the TE restoring force. This 
dramatic difference is due to the presence of electric field components that are perpendicular 
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to the dielectric boundaries. Specifically, at the left boundary of the output waveguide = 
plane), Ez induces a polarization surface charge density, while gives it a transverse kick. 
The result, Fx2,i, given by Eq. ([2]), is plotted in Fig. [2](b). This force is negative for the 
entire range of d and is seen to comprise almost all of the total force in the self alignment 
regime. Qualitatively, it may be associated with a dipole induced by the guided mode 
which has E^ oc sm{kxd), and the strength of the dipole per power increases with d as 
the mode confinement increases. At the z = interface, the dipole is roughly inverted, 
and consequently the two parts attract each other, in both the transverse and longitudinal 
directions. Moreover, this force grows rapidly as a function of A and is therefore responsible 
for the derivative approximation being less accurate than for the TE case. This is seen in 
Fig. El^b) where Frc2 is plotted for four different values of rf: d = 55 nm where instability 
is observed, d = 55 nm about where the derivative vanishes at A = 0, (i = 85 nm which 
exhibits self alignment, and d = 110 nm, where there is strong self alignment which is further 
discussed below. 

III. A WAVEGUIDE BROKEN BY AN OFFSET AND A GAP 

So far we have established the fundamental tendency of a waveguide to self-align or 
be in instability by considering zero gap and a perturbation in a continuous waveguide in 
the form of an offset. We next extend the discussion by introducing a longitudinal gap 
g ^ 0. Figure Hl^a) and Fig. ID^b) show contours of F^2 in the g-A plane obtained by FEM 
simulations, for TE and TM incident modes, respectively. The slab half-width is assumed 
to be c? = 110 nm, where according to Fig. [2l the TE mode places the system in instability, 
whereas the TM mode causes self alignment. In both cases the g = behavior extends 
to larger values of g, but the decay of the TE repulsive force with the offset and gap is 
much slower than that of the attractive TM force, as seen by the scales of g and A of the 
two frames of Fig. [2l The maximal TE force is about 0.4Fo and it is obtained for about 
A ~ 130 nm and (7 ~ 50 nm. A waveguide cantilever aX g = is, in fact, in a bistable 
state where a small offset may result in a deflection force that would eventually be balanced 
by the mechanical force. The attractive TM force is obtained for g = and A ~ 45 nm, 
and it is about O.TFq. We further found that a strong longitudinal force is pulling the two 
waveguides towards each other at a force of about 2Fq for (yf ~ 20 nm. For large enough 
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FIG. 4: Contours of Fx2/Fq as a function of g and A for d = 110 nm and = 3.48^. (a) TE 
incident mode, (b) TM incident mode. 

values of g the TM force becomes repulsive, corresponding to radiation pressure. 

The self alignment and the instability may be tested experimentally by fabricating on a 
chip two waveguide cantilevers with an offset and a gap. For instance, for Pin = 30 mW, 

= 100 pN, and at g = 20 nm we obtain Fx2 — O.I-Fq according to Fig. Mjo), which is 
about lOpN. This is of the order of magnitude of force that was shown to actuate a silicon 
cantilever [l^. The deflection of the waveguide in such a system may be viewed by the 
nonlinear input/output behavior, as the output power increases when the cantilevers tend 
to self align. Moreover, the mechanical effect is doubled by the fact that a similar force 
opposite in sign is exerted on both cantilevers. The two cantilevers may also be vibrated 
at their mechanical resonance by modulating the incident power, resulting in a system that 
may be suitable for applications such as sensing. 

IV. CONCLUSION 

In conclusion, we demonstrated a novel effect of light forces in the form of self alignment 
or instability of a waveguide broken by an offset and a gap. The waveguide size and mode 
polarization determine which of the two regimes the waveguide is in. Closed form expressions 
for the transverse forces were given for the case of a small offset and no gap. The forces 
described here are unique as they are due to the presence of both the guided mode and the 
scattered light from the discontinuity. Strong self alignment for a TM mode is caused by 
near field interaction of the polarization surface charges created by the longitudinal electric 
field. We are currently looking into the possibilities of an experimental realization that will 
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demonstrate the effects discussed here. 
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